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Abstract 



To measure an observable of a quantum mechanical system leaves it in one of 
its eigenstates and the result of the measurement is one of its eigenvalues. This 
process is shown to be a computational resource. It allows one, in principle, to 
■ diagonalize hermitean (N x N) matrices by quantum mechanical measurements 

only. To do so, one considers the given matrix as an observable of a single spin with 
appropriate length s which can be measured using a generalized Stern-Gerlach 
apparatus. Then, each run provides one eigenvalue of the observable. As it is 
based on the 'collapse of the wave function' associated with a measurement, the 
Q_i| procedure is neither a digital nor an analog calculation — it defines thus a new 

quantum mechanical method of computation. 

cr. 



Non-classical features of quantum mechanics such as Heisenberg's uncertainty relation 
and entanglement have intrigued physicists for several decades. From a classical point 
of view, quantum mechanics imposes constraints on the ways to talk about nature. An 
electron does not "have" position and momentum as does a billiard ball. Similarly, if 
a photon is entangled with a second one — possibly very far away — one cannot ascribe 
properties to it as is done for an individual classical particle. The lesson to be learned 
is that classical intuition about the macroscopic world simply does not extrapolate into 
the microscopic world. 

In recent years, an entirely different attitude towards quantum theory has been put 
forward. The focus is no longer on attempts to come to terms with its strange features 
but to capitalize on its both counter-intuitive and well-established properties. In this 
way, surprising methods have been uncovered to solve specific problems by means which 
have no classical equivalent: quantum cryptography, for example, allows one to establish 
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secure keys for secret transmission of information PJ; entanglement M is used as a 
tool to set up powerful quantum algorithms which do factor large integers much more 
efficiently than any classical algorithm ||. Throughout, these new techniques rely on 
the measurement of quantum mechanical observables as a reliable tool. This is also true 
for quantum error correction j|, [5| required to let any potential algorithm run. 

Here the purpose is to point out that the bare 'projection' effected by a quantum 
mechanical measurement does possess computational power itself. As will be shown 
below, it can be used to solve explicitly at least one specific computational task, namely 
to determine eigenstates and eigenvalues of hermitean (N x N) matrices. 

The diagonalization of hermitean matrices is a recurrent problem in mathematics, 
physics, and related fields. Using the notation of a quantum physicist the problem reads 
as follows. Given a self-adjoint operator A acting on a Hilbert space 7i of dimension N, 
one needs to determine its eigenstates \A n ), n = 1, . . . , N, and its A" real eigenvalues A n 
satisfying A \A n ) = A n \A n ), n = 1, . . . , N. If normalized to one, the eigenstates consti- 
tute a complete orthonormal basis of the space H: Y, n =i \ A n )(A n \ = 1, (A n \A n i) = 5 nn i. 
The standard solution from linear algebra H is to write down the eigenvalue equation 
with respect to a given orthonormal basis \k),k — 1, . . . , N, say. The A^ 2 matrix elements 
Afcfc' = (k\A\k') determine the operator A uniquely and its eigenstates are characterized 
by the coefficients (A n ) k = A nk in the expansion \A n ) = J2kA n k\k). The number A 
is an eigenvalue of A if the characteristic polynomial -Pa(A) of the matrix A vanishes, 
-Pa (A) = det (A — AE) = 0, where E is the (N x N) unit matrix. Once the A" roots A n 
of the polynomial Pa are known, the non-zero solutions of the equation 

(A-A n )A n = 0, n = l,...,N, (1) 

provide the eigenvectors \A n ) in the basis \k). Analytic expressions for the eigenvalues 
A n in terms of the elements of A exist only if A^ < 4. In general, numerical methods are 
required to determine approximately the roots of -Pa(A). 

The quantum diagonalization of hermitean matrices is based on the assumption that 
the behaviour of a spin s is described correctly by non-relativistic quantum mechanics. 
This method will make use of the 'collapse of the wave function' as computational 
resource. Note that the procedure does not depend on a particular interpretation of 
quantum mechanics. Five steps are necessary to achieve the diagonalization of a given 
matrix A (supposed for simplicity not to have degenerate eigenvalues). The individual 
steps will be described first in a condensed form; subsequently, commentaries explain 
the technical details. 

1. Standard form of A: Write the hermitean (N x N) matrix A as a combination of 
linearly independent hermitean multipole operators T„, v = 0, . . . , A^ 2 — 1, 

N 2 -l j 

A = J2 a - T - > ^=m Tt t AT ^ G R (2) 



2 



2. Identification of an observable: Interpret the matrix A as an observable hU for a 
single quantum spin S with quantum number s = (N — l)/2, 

N' 2 -l 

H A (S)= a*T„(S), (3) 

using the expression of the multipoles T V (S) in terms of the components of a spin. 

3. Setting up a measuring device: Construct an apparatus app(H A ) suitable to mea- 
sure the observable H^. 

4. Determination of the eigenvalues: Carry out measurements with the apparatus 
app(H 7 4) on a spin s prepared in a homogeneous mixture p = 1/ (2s + 1). The 
output of each measurement will be one of the eigenvalues A n of the matrix A. 
After sufficiently many repetitions, all eigenvalues will be known. 

5. Determination of the eigenstates: Calculate the eigenstates \A n ) of the matrix A 
on the basis of Eq. (|I]) and the experimentally determined eigenvalues A n . Al- 
ternatively, determine the eigenstates \A n ) experimentally by methods of state 
reconstruction. 

Thus, the matrix A has been diagonalized without calculating the zeroes of its char- 
acteristic polynomial by traditional means. The fourth step solves the hard part of the 
eigenvalue problem since it provides the eigenvalues A n of the matrix A. The comments 
to follow provide the background necessary to perform the individual steps. Emphasis 
will be both on the construction of a device measuring for a given hermitean operator 
(Step 3) and on the working of a quantum mechanical measurement (Step 4). 

Ad 1: The A^ 2 self-adjoint multipole operators T u = Tj, form a basis in the space of 
hermitean operators acting on an A^-dimensional Hilbert space 7i [/J. Two multipoles 
are orthogonal with respect to a scalar product defined as the trace of their product: 
(1/AO Tr [T U T U ,] = 5 W >. 

Consider now a Hilbert space 7i s of dimension (2s + 1) which carries an irreducible 
representation of the group SU{2) with the spin components (Si, S2, S3) as generators. 
Then, the multipoles T u ,u = 1,...,N 2 — 1, are given by the symmetrized products 
Sj x Sj 2 ■ ■ ■ Sj a ,ji = 1, 2, 3, and a = 0, 1, . . . , 2s, after subtracting off the trace (define 
To = T'W = E, the (N x N) unit matrix). The index a labels (2s + 1) classes with 
(2a+l) elements transforming among themselves under rotations; for the sake of brevity, 
a collective index v = (a;ji, . . . ,jk) is used. Explicitly, the lowest multipoles read 

T ( 1 ) = c. T (2) = ifs.S- +S- S- ) - ihllS- S (A) 

1 J -'31 3132 2^ 31 32 32 n> 3 31 32 ' ^ ' 

The set {T^} is a basis for the hermitean operators on 7i s . 
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Ad 2: Since the multipoles are expressed explicitly as a function of the spin compo- 
nents not exceeding the power 2s, it is justified to consider them and, a fortiori, the 
quantity as an observable for a spin s. 

Ad 3: It is natural to expect that every self-adjoint operator B comes along with an 
apparatus app(-B) capable of measuring it || . For particle systems, setting up such a 
device remains a challenging task for an experimenter. 

For spin systems, the situation is different, however. Swift and Wright |7[] have shown 
how to devise, in principle, a generalized Stern- Gerlach apparatus which measures any 
observable HU(S) — just as a traditional Stern-Gerlach apparatus measures the spin com- 
ponent n-S along the direction n. The construction requires that arbitrary static electric 
and magnetic fields, consistent with Maxwell's equations, can be created in the labora- 
tory. To construct an apparatus app(H^4) means to identify a spin Hamiltonian H(r,S) 
which splits an incoming beam of particles with spin s into subbeams corresponding to 
the eigenvalues A n . The most general Hamiltonian acting on the Hilbert space 7i of a 
spin s reads 

N 2 -l 

H(r,S)= J2 (5) 

with traceless (except for v = 0) symmetric expansion coefficients $„(r)(= 3>j*j 2 Jfe ( r )) 
which vary in space. Tune the electric and magnetic fields in such a way that the 
coefficients 3>i/(r) and its first derivative with respect to some spatial direction, r\, say, 
satisfy 

t , (r= „ ) = ^l = v (6) 

This is always possible with realistic fields satisfying Maxwell's equations. Then, the 
Hamiltonian in ([|) has two important properties, (i) At the origin, r = 0, it coincides 
with the matrix H^. (ii) Suppose that a beam of particles with spin s enters the 
generalized Stern-Gerlach apparatus app(hU) just described. At its center, particles 
in an eigenstate \A n ), say, will experience a force in the r\ direction given (up to second 
order in distance from the center) by 

w( ^ d(A n \H(r = 0,S)\A n ) 

Fx{r = 0) = = -A n , n = 1, . . . , 2s + 1 . (7) 

Consequently, particles with a spin projected onto one of the eigenstates \A n ) of the 
operator are separated spatially by this apparatus. The procedure is entirely analo- 
gous to that for a spin 1/2 where a familiar Stern-Gerlach apparatus is used (see for 
details). 

Ad 4: The 'projection postulate' of quantum mechanics describes the effect of mea- 
suring an observable Bona system S by means of an apparatus app(-B). If the system 
is prepared initially in a state with density matrix p one has: 

app(£?) : p (B n ;p n ) , p n = Tr [pp n ] . (8) 
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The action of the apparatus is, with probability p n , to throw the system into an eigen- 
state p n = \B n )(B n \ of the observable B\ the outcome of the measurement is given by 
the associated eigenvalue B n . By the way, the notion of 'collapse' or 'projection' can 
be avoided by characterizing the process indirectly by refering to "repeatable measure- 
ments" 0. 

The outcome of an individual measurement cannot be predicted due to the proba- 
bilistic character of quantum mechanics. Therefore, the probabilities p n , resulting from 
(infinitely often) repeated measurements on identically prepared systems, represent the 
essential link between theory and experiment. They provide information about the state 
of the system conditioned by the selected observable. Thus, a measurement reveals (or 
confirms) properties of the state p of the system while the observable B at hand is 
assumed to be known, including its eigenstates and eigenvalues. To put it differently, 
the observable defines the scope of the possible results of a measurement: the only pos- 
sible outcomes are its eigenvalues B n , and, directly after the measurement the system 
necessarily resides in the corresponding state \B n ). 

As the occurrence of the eigenvalues is purely probabilistic, one needs to repeat the 
experiment until all values A n have been obtained. If the spin s is prepared initially in a 
homogeneous mixture, p = E/(2s + 1), the (2s + 1) possible outcomes occur with equal 
probability. The probability not to have obtained one specific value A n after Nq ^> N 
measurements equals l/(2s + 1) , decreasing exponentially with Nq. 

Ad 5: It would be very convenient now to 'read out' directly the quantum state p n 
obtained from a single measurement with result A n . However, due to the no-cloning 



theorem ||T0| , pT| , an unknown state cannot be determined if only one copy of it is 
available. Upon repeating the measurement a large number of times and keeping only 
those states with the same eigenvalue A n , one produces an ensemble of systems prepared 
identically in the state p n . This is sufficient to reconstruct an unknown state since a 
density matrix p can be written as 

i N 2 

P=mT, P »Q^ N = 2s + l, (9) 

where the coefficient P a = (n M |/5|n^) is the probability to find the system in a coherent 
spin state |n M ). The operators Q^, p, — 1, . . . , N 2 , form a basis for hermitian operators, 
similar to but different from the multipoles T u [|12| . Thus, Eq. (|9]) parametrizes p by 
expectation values P M which can be measured by a standard Stern-Gerlach apparatus.- 

In sum, the basic ingredient of quantum diagonalization is the 'collapse' of the wave 
function projecting any state onto a randomly selected eigenstate of the measured ob- 
servable. Generalizations of this approach are expected to include the diagonalization 
of unitary matrices and the determination of roots of polynomials. 

Usually, a measurement is thought to confirm or reveal some information about the 
state of the system. Here, on the contrary, the idea is to learn something about the 
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measured observable instead. Why is this possible? It is fundamental to realize that the 
input required to actually measure A differs from the output of the experiment: for a 
measurement of A, the construction of an apparatus app(A) is sufficient which is possible 
withoutknowing eigenvalues and eigenstates of A. Necessarily, after a measurement par- 
tial information about the spectral properties of the observable A is available according 
to (|). This is due to the constraints (i) that the possible outcomes of measuring A 
are its eigenvalues and (ii) that the system subsequently will occupy the corresponding 
eigenstate. Thus, if the eigenstates and eigenvalues of A not known initially, one indeed 
acquires information about them by measuring A. 

The quantum mechanical diagonalization appears to be neither an analog nor a 
digital calculation. It is not based on the representation of a mathematical equation 
in terms of a physical system which then would 'simulate' it. Similarly, no 'software 
program' is executed which would implement an diagonalization algorithm. One might 
best describe the measuring device app(hU) as a 'special purpose machine' based on the 
projection postulate. 

For the time being, the method introduced here is important from a conceptual but 
not a technological point of view. On the one hand, the diagonalization of matrices is 
not a hard problem such as factorization of large integer numbers; on the other, the 
actual implementation in the laboratory is challenging. It is important, however, that 
there is no physical principle which would forbid the construction of such a machine. 
Further, it is expected to be fruitful from a conceptual point of view since it provides 
a different perspective on the projection postulate |13H . Quantum diagonalization as 
introduced here shows that — in an unexpected way — standard quantum mechanics at- 
tributes computational power to the measurement of an observable. The fact that one 
can use a measurement to perform calculations might turn into an argument in favor of 
the 'reality' of the quantum mechanical projection postulate. 
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